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ABSTRACT 

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedi- 
mentation of heavy elements. Motivated by observations of Type I X-ray bursts from sources with 
extremely low persistent accretion luminosities, Lx < 10 36 ergs s _1 (~ O.OlLEdd), we study how sed- 
imentation affects the distribution of isotopes and the ignition of H and He in the envelope of an 
accreting neutron star. For local mass accretion rates m < 10~ 2 mEdd (for which the ignition of H 
is unstable), where mEdd = 8.8 x 10 4 g cm -2 s _1 , the helium and CNO elements sediment out of 
the accreted fuel before reaching a temperature where H would ignite. Using one-zone calculations 
of the thermonuclear burning, we find a range of accretion rates for which the unstable H ignition 
does not trigger unstable He burning. This range depends on the emergent flux from reactions in the 
deep neutron star crust; for F = 0.1 MeV(m/m u ), the range is 3 x 10 _3 mEdd ^ ih < 10~ 2 77iEdd- We 
speculate that sources accreting in this range will build up a massive He layer that later produces an 
energetic and long X-ray burst. At mass accretion rates lower than this range, we find that the H flash 
leads to a strong mixed H/Hc flash. Surprisingly, even at accretion rates m > 0.1?riEdd, although the 
H and He do not completely segregate, the H abundance at the base of the accumulated layer is still 
reduced. While following the evolution of the X-ray burst is beyond the scope of this introductory 
paper, we note that the reduced proton-to-seed ratio favors the production of 12 C — an important 
ingredient for subsequent superbursts. 

Subject headings: diffusion — stars: neutron — X-rays: binaries — X-rays: bursts 



1. INTRODUCTION 

An ionized plasma in a gravitational field develops an 
electric field sufficient to levitate the ions and ensure 
overall charge neutrality. When there is more than one 
species of ion present, the ions will experience a differ- 
ential force: lighter ions float upward (defined by the lo- 
cal gravitational field) and heavier ions sink downward. 
The high surface gravity of compact objects makes the 
timescale for isotopes to stratify often comparable to or 
faster than other timescales of interest. In particular, 
sedimentation of heavier isotopes 4 is important in un- 
derstanding the surface compositions of whit e dwarf stars 
ijVauclair et alJ|1979UPaauette et al.ll986bl) and isolated 
neutron stars ijChang fe Bildstenll2003|) . For accreting 
white dwarfs, diffusion between the accreted envelope 
and underlying white dwarf was proposed as a means 
to enrich the ejecta of c l assical novae in CNO isotopes 
(IPriahrik fc KovetzllT98l llben. Fuiimoto. fe MacDonaldl 
11992ft . 

Accreting neutron stars, with their strong surface grav- 
ity w 2.0 x 10 14 cms -2 , are an ideal place to look for 
the effects of sedimentation. The sedimentation of heavy 
elements and resulting nucleosynthesis in the envelope 
of isolated neutron stars cooling from birth was first 
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described bv iRose nl (Il969 albf) and has b een studied in 
detail bv IChang fc Bildstenl {20031 [20041. For accret- 
ing neutron stars, the rapid stratification removes heavy 
nuclei from th e photosphere for accretion rates M < 
10 _12 M Q yr -1 ijBildsten. Sakoeter. fc Wassermanlll992y 
Deeper in the neutron star envelope, the differentiation 
of the isotopes can alter the nuclear burning, namely un- 
stable H/He burning and rapid proton-capture process 
(rp-process), that powers Type I X-ray bursts. Some 
estimates of the relative i mportance of sedimentation 
and diffusion were made bv lWallace. Wooslev. fc Weaver! 
l|1982ft . who studied accretion at M = 2 x 10 -11 M Q yr _1 
and argued that the partial separation of the H and He 
layers might play a role in setting the ignition conditions 
and subsequent burst nucleosynthesis. 

Recent long-term monitoring of the galactic cen- 
ter with BeppoSAX led to the discovery of nine 
"burst-only sources" (see ICornelisse et all 12004 
and references therein): SAX J1324.5-6313, 1RXS 
J1718.4-4029, GRS 1741.9-2853, SAX J1752. 3-3138, 
SAX J1743. 5-2349, SAX J1806.5-2215, 
SAX J1818. 7+1424, SAX J1828. 5-1037, and 
SAX J2224.9+5421. These sources did not have 
persistent fluxes detectable with the BeppoSAX /WFC, 
thus must have extremely low accretion luminosi- 
ties (Lx < lO^erg ss " 1 ^ 0.01L E dd)- Even earlier, 
iGotthelf fc Kulkarnil l)1997[) had detected a low- 
luminosity X-ray burst (peak burst luminosity w 0.02 
-^Edd) from the globular cluster M28. If the accretion 
is not concentrated onto a small surface area, so that 

t«@ftbvl«a c g^ is rh ~ M/(47ri? 2 ), then the 
sedimentation ilmescale for CNO nuclei, defined as the 
time required to move a scale height relative to the cen- 
ter of mass of a fluid element, is less than the accretion 
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timescale, defined as the time for a fluid element to reach 
a given depth. This is demonstrated in § where we 
show that the sedimentation timescale of CNO nuclei is 
less than the accretion timescale at m < 0.05 riiEdd- Here 
m E dd = 2m p c/[(l + X K )a TH R} = 8.8 x 10 4 g cm~ 2 s" 1 
is the local Eddington mass accretion rate for a solar 
composition, oth is the Thomson scattering cross 
section, and Ah is the hydrogen mass fraction. An 
estimate of t h is acc retion rate was noted earlier by 
I Wallace et al.l (jl982j) . As a result, sedimentation in 
the accreted envelope must be considered in treating 
the unstable ignition of hydrogen and helium for these 
low-m sources. 

Motivated by these observations, we explore in this 
paper the unstable ignition of hydrogen and helium 
at low accretion rates and pay particular attention 
to the regime where the mass accretion rate is less 
than the critical rate needed for stable H burning 
ijFuiimoto. Hanawa. fc Mivai3ll98H) . We study how the 
sedimentation of CNO nuclei affects the unstable ignition 
of H and He in an accreted neutron star envelope. At low 
accretion rates, for which the sedimentation timescale 
is comparable to or shorter than the burst recurrence 
timescale, the partial separation of the H and He lay- 
ers is important in setting the ignition conditions and 
subsequent burst nucleosynthesis. In particular, this af- 
fects the issue of whether unstable H burning triggers the 
X-ray burst. We find that there is a range of mass ac- 
cretion rates for which H burning is unstable (according 
to a simple one-zone prescription) but too weak to ig- 
nite He. We speculate that this may allow the formation 
of a deep, cold He layer, and we discuss how the out- 
come depends on the flux emergent from reactions in the 
deep crust of the neutron star. At accretion rates lower 
than this range, we find that the H flash leads to a strong 
H/He flash. Our work is directly applic able to the bursts 
obser ved from the "burst-only sources" l)Cornelisse et alj 
|2004[K for which the time between bursts is observed to 
be longer than the sedimentation timescale. 

At higher accretion rates, sedimentation is too slow to 
separate the accreted envelope; nevertheless, an abun- 
dance gradient develops. The concentration of H at the 
base of the accreted envelope is reduced, both because 
H tends to rise and CNO to sink and also because the 
H is consumed faster due to the enhanced CNO abun- 
dance. Thus, at the depth of burst ignition (triggered by 
the 3a reaction) the proton abundance is reduced. This 
may be important for determining the amount of 12 C 
produced, which is the most likely fu el for igniting su- 
perbursts l)Cumming fc Bildstenll200lD . We find that the 
elemental separation plays a role in setting the ignition 
conditions of the X-ray burst and the resulting nucleosyn- 
thesis even at higher accretion rates m > O.lmEdd- 

In this paper, we first provide, in § |2 a review of the 
formalism for calculating the elemental sedimentation. 
We then incorporate, in §|21 our calculations of the differ- 
ential isotopic velocities into a model of the accumulating 
neutron star envelope. With this model we examine (§EJ 
the relative abundances of H, He, and CNO isotopes at 
ignition, and we use a one-zone model to follow the sub- 
sequent thermonuclear burning. We conclude in § El by 
discussing directions for future research. 

2. SEDIMENTATION AND DIFFUSION 



To describe a multifluid gas, we follow the treatment of 
iBurgers! (^969) , which constructs the equations from suc- 
cessive moments of the Boltzmann equation. Our plane- 
parallel atmosphere is in hydrostatic equilibrium, and for 
simplicity we negle ct the terms coupling th ermal and par- 
ticle diffusion (see lPaauette et aTlll986ai and references 
therein). From these assumptions, we then have for each 
species s an equation of continuity and momentum con- 
servation, 



+ n s A s m u g - n s Z s eE = ^ K st (w t - w s ). (2) 

t 

Here species s has mass A s m u , m n being the atomic 
mass unit, charge Z s e 1 density n s , partial pressure P s , 
and velocity u s . We denote the substantial derivative for 
this species by D s /Dt = dt + u s d r . The force terms in 
equation J5J) result from the gravitational field g and the 
induced electric field E. The center-of-mass of a fluid el- 
ement moves with a velocity u = ^ s n x A K u s / y~l„ n^A*, 
and u = ml p to a good app roximation l)Bildsten et alJ 
ll99llBrown fc Bildstenll 19981) . The differential, or diffu- 
sion, velocity between species s and the fluid element is 
then w s = u s — u. These diffusion velocities satisfy mass 
and charge conservation, 

A s n s w s = 0, (3) 

s 

Z s n s w s = 0, (4) 

s 

with the electrons included in the summation. 

The right-hand side of equation J2J is the collision term 
between species s and t, with the resistance coefficient 
K st — n s n t {a s tv s t) being the velocity-weighted cross- 
section and v s t the center-of-mass relative velocity be- 
tween particles of types s and t. The nature of the ionic 
interaction is characterized by 

(ZV _{Z^_ ( p \ 1/3 / 10 8 K \ 

ak B T (A) 1 / 3 Vl0 5 gcm-V V T J ' 

(5) 

where a is the average ion spacing, defined by 
4?ra 3 n/3 = 1 with n = pN A /(A). The ions 
are not strongly interacting, and the temperature 
is sufficiently high that the ion quantum occupation 
is small, n/n Q = n[2irh 2 /{Am u k B T)} 3 ^ = 3.2 x 
10-^-5/2(^/105 g C nr- 3 )(T/10 8 K) 3 / 2 < 1. The 
ions are not as s trongly interacting as in cooling 
white dwarfs (see IDelove fe Bildsterjr i2"002. and ref- 
erences therein). The unstable ignition of H and 
He occurs where T < 1, which is where the resis- 
tance coefficients K are most uncertain. Appendix 
lAl describes in detail our choice for computing the 

resistance coefficient s; in short, we use the fit of 

iFon^aine fc Michauc)| l|1 979[). which agree s with other fits 
jMuchmorel 11984 iPaauette et alJll986aft in the regime 
r ~ 1 , has the correct scaling at V ^> 1 IIBildsten fc Hall 
120011 IHansen. Jolv. fc M cDonaldl I1985|) . and goes, in 
the limit T <C 1, to K® t (see Appendix computed 
using a Coulomb po tential with a long -range cutoff 
l)Chapman fc Cowlindli"952l lBurgersjll969() . 
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Before describing our detailed calculations of the ac- 
creted envelope structure, we present a simpler case to 
illustrate how sedimentation changes the structure of the 
envelope. In the case of a trace species, labeled 2, in a 
background of a species labeled 1, i.e., n\ 3> n-i and 
w\ = 0, equations (J2J separate: the right-hand side of 
the equations vanishes, thereby fixing the electric field, 
eE = m u gAi/(Zi + 1) where the electrons are non- 
degenerate and eE = m n gA\/Zi where the electrons are 
degenerate. Substituting E into the equation of motion 
for species 2 then determines the sedimentation velocity, 
ii; sc d = w 2 — n 2 (A 2 m u g — Z 2 eE)/Ki 2 . We chose the sign 
of w se d to be positive if species 2 moves downward. 

Using the Stokes-Einstein relation to determine 
K-\2 from the d rag c oefficient for a liquid sphere 
(Bildsten fc Halll 120011: see Appendix EJ, and a 
nonrclativistic electron equation of state, we find 
the sedimentation veloci ty of a trace nucleus (see 
iBrown. Bildsten. fc ChaTiell200l 



w sed = 2 x 10" 



Tf 3 (A 2 Z X 



5l4- .6 
f5 



Z 2 A 1 )A° 1 1 



where we 

-i!4 



(6) 



cms 

-3 



use 

-2 



the common shorthand <?i4 = 
), T 7 = T/(10 7 K), and p$ = 
In a pure H plasma, the sedimen- 



sAio 1 

p/(10 5 gcm- 3 ) 

tation velocity 5 is greater than the mean velocity 
u = rh/p for m < 400 g cm" 2 s" 1 ^ 3 ^- 4 ^ - Z 2 )/Z$- 3 ; 
under conditions at which H ignites (T7 ~ 5; ps ~ 4) 
this corresponds to m < 0.02rriEdd and m < 0.05mEdd 
for a trace 4 He and 12 C nucleus, respectively. 

For a more refined estimate, we transform to a frame 
co-moving with a fluid element and solve for the displace- 
ment of a test particle from the center-of-mass of the fluid 
element, 



dL 



sod 



dL 



dy 



scd . — l 

m 



df 



w scd m 



(7) 



while simultaneously solving for the thermal and com- 
positional structure of the accreting envelope (see § |3| 
eq. |5] and ^U])- Here dy — —pdz — rhdt is the col- 
umn density. The electric field is compu ted using an 
analytical formula ijChang fc Bildsten 2003) that is valid 
for arbitrary electron degeneracy. Figure ^ shows L so d 
{thick lines) for a trace 4 He nucleus in a bath of H as 
a function of the column y. For comparison, we also 
plot the pressure scale height Hp = P/(pg) (thin lines). 
Both L se d and Hp are shown for 4 different values of 
rh: 100 (solid lines), 500 (dotted lines), 2500 (dashed 
lines) and 10 4 g cm~ 2 s _1 (dot-dashed line). For later 
convenience, we define y C ross(m) to be the column den- 
sity at which L se d = Hp. An accreted fluid element 
will be significantly stratified by the time it reaches a 
column Across- For a fluid mixture, the growing con- 
centration gradient slows the drift velocity as diffusion 
becomes important. The balance of diffusion and sedi- 
mentation sets the equilibrium abundance profile, and so 
over the time to reach a column y croS s the fluid element 
roughly comes into diffusive equilibrium. y cr oss(fn) is a 
strong function of accretion rate and surface gravity: for 
a solar composition plasma and the velocity estimate for 

5 Throughout this paper we use a Newtonian metric and assume 
a neutron star mass 1.4 Mq and radius 10 km, so that 314 = 1.9. 



J3 



c 

s 

'■3 

C/3 



10 4 



10 2 r 



10' 



10" 



10- 



10- 







r 

m = 100 . 

- - - ' "- 


_ . • ■ ' — 

• ■ ' ' — :: ^ ; . 


: . „- 




/ ...-•'500^^ " ' 








If 














[1/ 




if 




ii 1 1 1 





10 3 10 4 10 5 10 6 10 7 
Column Density [g cm 2 ] 



MY 



Fig. 1. — Sedimentation depth L so( j for 4 He (thick lines) and 
pressure scale height H p (thin lines) as a function of column density 
for 4 mass accretion rates: rh = 100 (solid lines), 500 (dotted lines), 
2500 (dashed lines) and 10 4 g cm -2 s -1 (dot-dashed lines). 



a trace He nucleus in H (eq. 0) and the scale height 
for a degenerate, nonrelativistic electron gas, y croS s ~ 
2.4 x 10 8 m3 5/6 T 7 ~ 5/4 ffl4 31/6 gem" 2 . For trace CNO ele- 
ments in H, y cross ~ 10 7 rh 2 ^ 7 ' 6 1 4 gT^' ' 6 gem -2 . Note 
that with the approximations in eq. ©, the sedimenta- 
tion velocity of trace CNO elements in He is zero and we 
must use the full equation of state, including Coulomb 
effects, to compute w so d- 

For HCNO burning, H is consumed on a timescale « 
5 x 10 4 (0.02/Z C no)(^h/0.71) s, or equivalently at a 
column 



2/bm„ = ~m ( ^ ) = 3x10 s gem" 



0.1 m E dd 



0.02 



,ZcNO 

(8) 

Using the previously defined y C ross, we determine the 
mass accretion rate below which the CNO elements are 
separated from H by setting y cross = y^um- We find the 
critical accretion rate at which the CNO nuclei and H 
nuclei separate on a timescale comparable to them to 

ignite, m C rit/TOEdd « 0.02T 7 15 ^ 38 (^y 19 . For typical con- 
ditions for H ignition, T7 ks 5, and using 514 = 1.9, we es- 
timate that CNO elements are depleted before H ignites 
for rh < O.lmEdd- Helium sinks more slowly than CNO; 
in the absence of CNO sedimentation, He would settle 
out of the H layer prior to H ignition for rh < 0.06 mEdd- 
These simple estimates will be a useful guide for our more 
detailed solutions, as we now describe. 

3. EVOLUTION OF AN ACCRETING ENVELOPE WITH 
SEDIMENTATION 

We now build on the simple estimates of section [21 and 
numerically solve for evolution of an accreting envelope 
prior to unstable ignition. In addition to equations l(T)l 
and J2J we have the equations for temperature and heat 
flux, 



dT _ F 
dy pK 

— -C (— rh— 

dy P \dt dy 



CpTrh 

y 



(9) 



V ad - £nucl ,(10) 
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where Cp is the specific heat at constant pressure and 
V a d = (d hxT/d In P)s- Note that we have neglected 
heat transport carried by the rising and sinking isotopes, 
as heat transport by photons and semi-degenerate elec- 
trons will be much more efficient. The compressional 
flux (terms proportional to m in eq. II O p are of order 
~ keV m/m„, at these accretion rates (jBrown fc Bildstenl 
f 998) and are smaller than the flux from steady H burn- 
ing (for m > O-Olrh-Edd) or from electr on captures in the 
crust (> 0.1 MeVm/m„ : lBrownll2000j) . 

We use a radiative-zero outer boundary, F = 
4/(3«y) [ctrT 4 - F acc ] , where F acc = GMm/(2R) is 
the accretion flux and n(p[y, T], T) is the opacity. For 
the inner boundary, at the bottom of the ash layer, 
we set the flux to a constant value. This flux is set 
by reactions in the deep crust — electron cap tures, neu- 
tron emissions, and pycnonuc lear reactions (Satol lT979T: 
lHaensel fc Zdunikl 11990. 12003]) — and for w > 0.1 m Ed d 
is roughly F = 0.1 MeV m/m u l|Brownll200f)(k at lower 
accretion rates the heat per nucleon flowing out is in- 
creased. Here we have computed solutions for different 
values of Fm u /m (see §0J). 

To simplify our calculations, we use the dispar- 
ity between the thermal and accretion timescales, 
Cpy 2 / (pK) <C y/m. At each timestep we solve equa- 
tions © and fTTijl with dT/dt — > 0, we presume that the 
envelope evolves through a sequence of steady-state ther- 
mal pro files. This technique has been used in previous 
studies l|Cumming fc Bildstenl l2000l iCooper et alll2005|) 
of H/Hc ignition. 

3.1. Microphysics of the Fuel Layer: Equation of State, 
Thermal Conductivity, and Nuclear Reactions 

Closing equations ©, and (fTUft requires a 

specification of the equation of state, the thermal con- 
ductivity, and the nuclear heating rate e nU ci- For 
the t hermal prope r ties, we mostly fol low previous pa- 
pers l|Brownl 12000: Br own et all 120021) . and here shall 
just review our choic es. We use a tabulatio n of the 
Helmholtz free energy (|Timmes fc Swestvll2000|) to com- 
pute the electronic EOS. We compute the ionic free en- 
ergy for the liquid phase, 1 < T < 175, from the fit of 
iChabrie r fc Potekhin] (|1 998f) and for the solid phase from 
the fit of lFarouki fc Hama guchl (|1993f) . The free ener- 
gies of these phases are equal at T = 178. We assume 
throughout that the ions are classical, i.e., we ignore 
quantization of phonon modes. For the temperatures 
and densities of interest, this is a good approximation. 

The relevant opacities for the temperatures and densi- 
ties in the neutron star envelope are Thomson scattering 
and free-fre e absorption. For free-free opacity, we use 
the fit from lSchatz et all l)1999f) which is reasonably ac- 
curate (fractional e rrors ~ 10 %) when compared agains t 
the calculations of |Itcfli 1 Jiaka£awa J _fc_I ^hvamal (l985); 
lltoh. Kuwashima. Ichihashi. fc Mutohl (|1991|t . We cal- 
cula te the Thomson scat tering opacity by using a 
fit ijBuchler fc Yuehl 119761) that reproduces the non- 
degenerate limit USampsonl IT951 and includes correc- 
tions for the relativistic and degenerate electrons. For the 
electron thermal conductivity, we write the relaxation 



(1980). We write the relaxation time for electron- ion 
scattering as 



time as r 



Here the term for electron- 
electron scattering, T s . n is compute d from the fit of 
iPotekhin et alJ l)1997j) to the results of lUrpin fc Yakovlevl 



47re 4 



pf,v F {A)m u 



((Z) 2 A e 



QA e ,imp) ■ (11) 



Here A ej i on is the Coulomb logarithm for a single ion 
of charge number (Z) and mass number (^4); Q = 
(Z 2 ) — (Z) 2 measures the impurity concentration; and 
A e ,imp is the Coulomb logarithm computed under the 
assumption that the ions are randomly distributed 
(lltoh fc Kohvamal 119931). For A „ ion we use the fitting 
formula of IPotekhin et al. (1993). 

To follow the nuclear burning, we use a reaction net- 
work with 22 isotopes: 1_2 H, 3 " 4 He, 7 Li, 7 Be, 8 B, 
12-13 Cj is-16 Nj u-i8 0| 17-19F, and 18 " 19 Ne. The light 
isotopes are included to cover hydrogen burning through 
p-p chains that are important at the low temperature 
end (m < 5.7 x 10" 4 m E dd = 50 g cm" 2 s" 1 ). The re- 
actio n rates are taken from t h e compilation REACLIB 
(see iThielemann et al.l 119861 iRauscher fc Thielemann! 
120001 and references therein), and most of them are 
origin ally from the compilation of iCauehlan fc Fowler] 
(1988). These rates are formally fitted for T > 10 8 K, 
but have the correct asymptotic form for the low- 
temperature reactions of interest here. We neglect 
electron screening. The rates for which experimental 
data were not avail able are typically taken from shell- 
model calculations (Rauschcr & Thielemann 2000). We 
include weak rates from iFuller et alJ l)1982f) and from 
ILanganke fc Martinez-Pinedol l|2001f) : weak interactions 
are not important, however, under conditions of interest 
in this paper. 

3.2. Numerical Solution for the Diffusive Velocities 

For a mixture of N ion species there are N + 
2 unknowns, the N + 1 velocities w s of the ions 
and electrons and the electric field E; for these un- 
knowns there are N + 2 equations (eqs. 0-01) ■ In- 
spection of equation J5J) reveals the velocities w s to 
contain an advective part depending on g, E, and 
d r T, and a diffusive part depending on d r n s . The 
mixed advective-diffusive nature of the equations makes 
simple centered-difference schemes, such as Crank- 
Nicho lson, unstable (see, for example, iFerziger fc Perioi 
2002). For num erical stability, we adopt a technique 
implemented by llben fc MacDonaldl l|1985j) and used 
in studies of the hydrogen-shell flash in evolved stars 
ijAlthaus. Serenelli. fc Benvenutoll200lD and pulsational 
instabilities in white dwarf s tars ijGautschv fc A lthaus 
l2002tlAlthaus fc C6rsicol2004() . The ion and electron dif- 
fusion velocities and the electric field are separated into 
two pieces, 



w t = to? 



d In nj 
dr 



E = EB ~ E Pi^r- 



(12) 
(13) 



j— ions 



where w 9 and E 9 are the components due to the grav- 
itational and electrical forces, and the sums on the 
right-hand side are only over the ions. Throughout 
this paper, we ignore the effect of thermal diffusion, as 
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this term is negligible unde r the conditions o f interest 
llStevenson fc Salr>eteHlT977t IMuchmore l 11 984ft . Insert- 
ing equations lfP2l and (|T~S|> into equations we 
use the fact that these equations must hold for arbitrary 
number density gradients for each species to equate the 
coefficients of dlnnj/dr\j—x t N- This gives us N(N + 2) 
equations for the variables (cr,j,/3j), in addition to the 
equations for the N + 2 variables (w 9 s ,E 9 ). We solve 
these linear equations by direct matrix inversion. 

With the <7y and [3j determined, we insert equa- 
tion H12|) into the isotopic continuity equation (eq. p^) 
and obtain 

— = ^(pw i Y i ) 



Dt 



d_ 

dy 



pYi 



Ed In rij 



(14) 



where Yi = riif (pN&) an d U = / pdz is the column den- 
sity. We solve for the Yi in each spat ial zone with a semi- 
impli cit finite difference method ijlben &: MacDonaldl 
1985). We tested this scheme by defining a closed box 
with an initially uniform mixture of H/Hc and an isother- 
mal temperature profile. We evolved the box in time 
until it reached diffusive equilibrium, and the resulting 
chemical pro file was in goo d agreement with the analyt- 
ical result of lAlcockl (^80, appendix). 6 

3.3. Evolution to Ignition 

We now implement the Lagrangian expression for the 
abundance evolution (eq. |14|) in a model of the accu- 
mulating neutron star envelope. We take the compo- 
sition of matter accreted from the companion star to 
be roughly solar composition in 1 K and 4 He, and dis- 
tribute the remaining mass evenly between 12 C and 16 0: 
X^H) = 0.71, X( 4 He) = 0.27, X( 12 C) = 0.01, and 
X( 16 0) = 0.01. Below the fuel layer are ashes from pre- 
vious X-ray bursts; we set the ash composition to 64 Zn, 
consistent with the findings of recent o ne-dimensional 
calculations of repeated X-ray bursts fWooslev et alJ 
l2004HFisker et all2005a|) that find A w 60 to be favored. 
We do not follow mixing between the fuel and ash layers 
and therefore this simple mo del neglects effect s arising 
from "compositional inertia" ijTaam et a!Jll993|) . 

For simplicity our scheme employs a fixed number of 
Lagrangian zones, so that our grid follows the advected 
fluid elements. A computation domain is closed, i.e., we 
set the diffusion velocities to zero at the boundaries. Our 
scheme is not as refined as schemes that add zones to the 
computational domain. The error added to this scheme is 
that H tends to accumulate at the upper boundary, which 
inhibits further sedimentation. Our results therefore err 
in the direction of a lesser separation of H and CNO, 
which is adequate for our initial exploration. 

We determine the size of a timestep by choosing the 
minimum of the thermal timescale, C p y 2 /(pK), the 
sedimentation timescale, Ay/(pw s ), and the burning 
timescale (|Y/F| nuc ). At each timestep we advect the 
matter inward and increase the pressure in each zone. 
We then update the composition using the reaction net- 
work followed by a diffusion step (eq. C21 _ EI)- Finally, 



6 N ote that the terms exp(— 5x/3) in equation (A14) of lAlcocU 
( 1980) should be exp(-5a;/2). 



we update the thermal profile by solving equations 
and (|10|l via relaxation l)Press et alJll992li . This process 
is repeated at each timestep until a thermal instability 
develops. We determine this ignition point using a simple 
one-zone criterion, (de nuc i/dT)p > (de coo i/dT)p, where 
£nuci is the nuclear heating rate and e CO oi = pKT/y 2 
is an approximation of the local cooling rate. Here K 
is the thermal conductivity. This scheme quantitatively 
agrees with linear stability analyses on ignition column 
den sities at regions where both methods find instabil- 
ity ijNaravan fc Hevll 12001 iCooner et"aT] l2005|h and it 
also roughly agrees wit h the ignition conditions found by 
multiz one calculations iWooslev et alJ2004t iFisker et alJ 
I2005a|> . 

4. RESULTS 

Using our accumulation model, we explore the evolu- 
tion of the accreting neutron star envelope for accretion 
rates 10 _4 -0.6 mEdd (10 g cm -2 s _1 -5 x 10 4 g cm~ 2 s _1 
for accretion of a solar composition mixture). For each 
accretion rate, we ran calculations both with and with- 
out isotopic sedimentation, and we follow the thermal 
and chemical evolution of the mixture until the envelope 
becomes thermally unstable, as described in §|3J W e firs t 
describe the qualitative features of our results (§ 14. ip . 
We then apply our findings to two problems: the phe- 
nomenon of X-ray bursts at low accretion rates (§ I4.2|l . 
and the production of 12 C to fuel subsequent superbursts, 
appropriate for sources accreting at rh > 0.1 mEdd (§E3l- 

4.1. Properties of the Accreted Envelope 

At accretion rates for which the temperature is hot 
enough for hydrogen to burn stably, hydrogen is removed 
from the base of the accreted fuel both by floating up- 
wards and by capturing onto the enhanced abundance of 
CNO nuclei that have settled downward. Figure|5]shows 
the mass fraction of H (solid lines), He (dotted lines) and 
CNO nuclei (dashed lines), relative to its initial value, 
at the base of the newly accreted layer for a local mass 
accretion rate 0.11 rfiEdd (10 4 g cm~ 2 s -1 ). Hydrogen 
ignition occurs after 10 3 s. The mass fraction is plot- 
ted against the time since arrival of the fluid element at 
the photosphere, so that the column corresponding to 
a given time t is y(t) — mt. We show for comparison 
solutions with diffusion (thick lines) and without diffu- 
sion (thin lines). At this accretion rate, H is consumed 
via the HCNO cycle (Fig. [2J , which increases the abun- 
dance of He. When sedimentation is included, the CNO 
mass fraction has doubled at the point where the H is 
consumed by HCNO burning. As a result, hydrogen is 
consumed on a shorter timescale, releasing a greater heat 
flux and increasing the He abundance. At the end of the 
calculation, the triple alpha reaction has started, increas- 
ing the number of CNO seed nuclei. We terminate the 
curves where the triple alpha reaction becomes thermally 
unstable; as is clear from the figure, the ignition column 
decreases by w 20% when sedimentation is included. 

This effect is even more pronounced at the accre- 
tion rate for which H docs not burn stably. Figure 01 
shows the same calculations for a lower accretion rate 
rh = 5.7 x 10~ 3 m Edd (500 g cm -2 s _1 ). For this calcu- 
lation, the base of the accreted envelope reaches diffusive 
equilibrium in 8 x 10 3 s (where the abundance curves flat- 
ten). This agrees with our simple estimates from § El (sec 
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Fig. 2. — The evolution of mass fractions of H (solid lines), 
He (dotted lines) and CNO elements (dashed lines), normalized to 
their initial values as a function of Lagrangian time t = y/rh for 
rh = 0.11 riiEdd- We show cases when diffusion is (thick lines) and 
is not (thin lines) included. The curves terminate where He ignites 
unstably. 




Time [s] 

Fig. 3. — Changes in mass fractions at the base of the fuel layer, 



for in = 5.7 x 10" 
those in Fig. [2] 



TOEdd- The curves have the same meaning as 



Fig. [3 and following discussion), namely the time for the 
envelope to be stratified is w Across/"" 1 ~ 7000 s. Note 
that as the envelope becomes stratified, the timescale for 
CNO to sink becomes set by its drift time in a He-rich 
plasma, which is considerably longer than the timescale 
for a H-rich plasma (eq. [5]) and becomes comparable to 
the timescale for He to sediment. This is why the curves 
in Fig. |21 all change slope at roughly the same time. 

We can compare the asymptotic values of the abun- 
dances in Figure wi th t h e ana lytical diffusive equi- 
librium solution of lAlcockl l)1980t appendix) . We use 
an ideal gas equation of state for simplicity. For a 
H/He plasma with / X^R) dz = 0.7 and / X( 4 Re) dz = 
0.3, the equilibrium solution at the base z = has 
X^R) = 0.29 and AT( 4 He) = 0.71, ie., the H has de- 
creased to 0.4 of its original abundance. For 12 C in a 
H-rich plasma with mass fraction of X( 1 H) = 0.98 and 
X( 12 C) = 0.02, we get X( 12 C,z = 0) = 0.23. If we 
now approximate the He-rich layer by a box containing 
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Fig. 4. — The temperature evolution of the base of the accreted 
layer as it is advected to deeper column. The tracks (solid lines) 
correspond to different local mass accretion rates and are in units of 
gem - s — 1 . The ignition curves of H and He when sedimentation 
is (dashed lines) and is not (dot-dashed lines) taken into account 
are shown as well. 



He and C with J X( 12 C)dz = 0.07 (so that the ratio of 
C/He matches 0.02/0.28), the equilibrium solution has 
X( 12 C, z = 0) = 0.17; that is, sedimentation increased 
the abundance of 12 C by a factor 2.4. Our numerical 
solution has X( 12 C, z = 0) increased by a factor w 4.5, 
which is somewhat larger than the analytical solution for 
12 C in 4 He, but much less than that for 12 C in H. 

The temperature evolution of our mass shell at the 
base of the accreted envelope is shown in Figure 0] Each 
curve is labeled by its local mass accretion rate in units 
of g cm" 2 s _1 . The x-axis indicates the column density; 
the temperature at the base of the accreted envelope fol- 
lows the trajectory in (P = yg,T) space shown in the 
figure. In all of the calculations shown, we include sedi- 
mentation. As with Figures|21and|31 the curves terminate 
where the envelope becomes unstable to a thin-shell in- 
stability (estimated from a one-zone calculation). The 
locus of these ignition points are indicated by the curves 
labeled "He ignition" (for rh > 2 x 10 3 g cm -2 s _1 ). and 
"H ignition" . Our conditions for unstable He ignition 
agree (in the absence of sedimentation and using an emer- 
gent flux from the crus t F = 0.15 MeV(m/m„)) ro ughly 
with those obtained bv lCumming fc Bildstenl l)2000|) . At 
rh = 0.1 riiEdd, our ignition column is 30% larger, but 
our temperature, hydrogen, and helium mass fractions 
agree to within 7%. This difference is likely caused by 
the H b urning rate being slower tha n the HCNO limit 
used bv ICumming fc Bildstenl pOOOj) at T 8 < 1.7. At 
these temperatures, the reaction 13 N(p, 7) 14 does not 
entirely dominate over the /3-decay branch, so that the 
rate is not entirely set by the decays of 14 and 15 0. The 
longer /3-decay time of 13 N (10 min halflife) decreases the 
total rate of H burning from the HCNO limit. 
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TABLE 1 
Critical Mass Accretion Rates 



Reference 


m c 2 
(g cm -2 s _1 ) 


(g cm -2 s — 1 ) 


this work (no sed.) 
Hanawa & Fuiimoto (1982) 
Bildsten (19981 
Naravan & Hevl (20031 
this work (with sed.) 


10 3 

3.2 X 10 2 

1.3 x 10 3 
3 x 10 2 
10 3 


2 X 10 3 

3 x 10 3 
8.6 x 10 3 
3 x 10 3 

5 x 10 3 



Note. — This comparison is for M = 1.4 Mq, 
R = 10 km, X = 7, and Z = 0.02. The re- 
sult of Naravan fc Hevll <2?X)3Tl is taken from the case 
of core temperature 1 s K and R = 10.4 km. For 
INarava n & Hevl (200i§), we interpret m c 2 as being the 
critical mass accretion rate required for prompt hydro- 
gen bursts. 



As noted by iFuiimoto et al.1 (0,981), there are three 
regim es of burning parameterized by m (see iBildstenl 
119981 and references therein): 1) tin < m C 2, for which 
hydrogen burns unstably; 2) m C 2 < rh < m c i, for 
which hydrogen burns stably and is completely consumed 
prior to unstable He ignition; and 3) m > m c i, for 
which hydrogen burns stably and is only partially con- 
sumed prior to unstable He ignition. In the absence 
of sedimentation, we find m C 2 ~ 10 3 g cm" 2 s , and 
m c i w 2 x 10 3 g cm -2 s _1 (see Tableland Fig.0. 

When sedimentation is included, m C 2 is unchanged:; 
the abundance of H at the base of the accreted envelope 
is depressed, however, for rh < m c i (see Fig. El - More- 
over, for accretion rates m < m c \ = 5 x 10 3 g cm -2 s _1 , 
He ignites in the absence of H. We emphasize, however, 
that the temperature at ignition and the total mass of H 
in the envelope is only slightly affected by sedimentation. 
The characteristics of the burst will depend on the inter- 
play between the t hermal instability and the growth of 
the co nvective zone ijWooslev et al.l2004llWeinberg et all 
2006); such a study is beyond the scope of this paper, 
but is clearly a crucial future step for understanding the 



burst physics. It is tantalizing that the accretion rate 
m c i at which mixed H/Hc ignition occurs is increased by 
a factor of 2 when sedimentation is taken into account, 
and we speculate that this might alleviate the discrep- 
ancy between the pre dicted transition in burst duration 
iFuiimoto et al.lll981fl and rece nt observations (see, for 
example Iden Hartog et al.ll2003j) . 

4.2. Bursts at Low Accretion Rates 
4.2.1. Observations 

As discussed in § ^ X-ray bursts with ex- 
tremely low persistent luminositie s ( < 1 36 ergs s ~ 1 ) 
have been discovered recently (see iCocchi et aljll99fi . 
200lt iKaptein et alJ l2000l iCornelisse et alJ [2002b: 
Arefiev fe AleksandrovicU l2004|) . In Table H we list 
the burst duration and persistent luminosity of all such 
known burst sources. Several of these did not have persis- 
tent fluxes detectable with the BeppoSAX /WFC and are 
known as "burst-only sources" . Follow-up observations 
with Chandra revealed that the sources' persistent lumi - 
nosities are 10 32 -10 33 ergss -1 l|Cornelisse et all 12002a). 
which is consistent with these sources being X-ray tran- 
sients. These bursts are very rare: on average, there is 
only one burst detected for every 10 6 s seconds of Bep- 
poSAX /WFC observation time. The bursts have rather 
short durations, < 20 s. At somewhat higher accretion 
rates, rh ~ O.OlmEdd, several bursts have been observed 
with longer durations, < 1000 s, that are intermediate be- 
tween mixed H/He bursts and superbursts (see the lower 
part of Table©. 

We note for completeness that still long er bursts were 
exhibi ted by sources at L < O.lLEdd- iCzernv et alJ 
l|1987l) observed (with Einstein /MPC) a burst of du- 
ration w 2500 s, when the pers istent lumino s ity w as 
2.5 x 10 37 eiRss- Hd/5 kpc) 2 (|Czernv et all Il987|) 7 . 
Fushi kTet al.1 l)1992l ) interpreted the burst as due to elec- 
tron captures at great depth. This burst had an initial 
e-folding time ~ 100 s, and was followed by a long tail. 
In morphology its lightcurve resembles that of a re cent 
burst observed from GX 3+1 l)Chenevez et al.ll2005|L 
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Fig. 5. — Mass fraction of hydrogen at the column where either 
H (open symbols, solid lines) or He (filled symbols, dotted lines) 
unstably ignites, as a function of mass accretion rates. We show 
results for which sedimentation is ignored (circles), and for which 
it is included (squares). 



4.2.2. Weak Hydrogen Flashes 

When the temperature in the neutron star envelope is 
sufficiently low, the CNO cycle becomes temperature de- 
pendent; as a result, the ignition of H becomes thermall y 
unstable at low accretion rates l)Fuiimoto et al.l ll9811. 
As a first investigation of the unstable ignition of H when 
the atmosphere is stratified, we perform several calcula- 
tions of the burst nucleosynthesis. We approximate the 
cooling by a one-zone finite differencing over the enve- 
lope, 

dT 



dt 



^-nucl ^-cool; 



(15) 



where e coo i = pKTy and we 
a reaction network. Included in this network are 
686 isotopes coveri ng proton-rich nuclei up to Xe (see 
ISchatz et all [2001) . The reaction rates are taken 
from the compilation reaclib (see Thielemann et al. 



1986; R auscher fc Thielemann! 120001 and references 



7 We adopt the fi ducial distance 5 
IRutledee et~al] 12001b]) 



kpc recommended by 



8 Peng et al. 



TABLE 2 

Bursts from Sources with Low Accretion Rates 



Source 


Burst Duration 


Persistent luminosity 


Reference 




c-folding time (s) 


(10 36 ergs -1 ) 




SAX 1324.5-6313 


6.0 ±0.1 


< 4 x 10~ 4 


C02, A04 


1RXS J1718.4-4029 


47.5 


0.07 


K00, A04, 105 


GRS 1741.9-2853 


8.8,11.0,16.0 


0.2 


C99, A04 


SAX 1752.3-3138 


21.9 


< 3 x 10" 4 


C01, A04 


SAX 1753.5-2349 


8.9 


< 4 x 10~ 4 


198, C02, A04 


SAX J1806. 5-2215 


4.0,9.0 


3 


198, A04 


SAX 1818.7+1424 


4.5 ±0.1 


< 4 x 10~ 4 


C02 


SAX 1828.5-1037 


11.2 ±0.6 


0.034 


C02, A04 


SAX J2224.9+5421 


2.6 ±0.2 


< 1.2 x 10" 3 


C02 


1RXS J1708. 5-3219 


~ 300 


1.5 


105 


SLX 1737-282 


~ 600 


< (0.5 - 1.8) 


102 


SLX 1735-269 


~ 700 


~ 4 


M05 



Referenc es. — A04:IAreiiev He AleksandrovichJ 120041) : C99, COEICocchi et all 
il99Sll200^: C02: ICornelisse et alJ 12002bl) ; K 00: IKaptein et alj 120001) : 198. l6T 
105: lin 't Zand et alJ I199SL 120021 12005al) : M05: IMolkov et all 12003) . 



therein), and consist of experimental rates and Hauser- 
Feshbach calculations with t he code non-smoker 
l)Rauscher fc Thielemann! 12000) . The initial tempera- 
ture, column density, and composition are taken from 
the values of the quasi-static calculation at the ignition 
point. For the unstable ignition of H, the temperature 
in the accreted envelope is set by the flux emergent from 
the deeper ocean and crust. As a result, the accretion 
rate at which the burst behavior changes will depend on 
assumptions about the heating in the crust. 

We shall first describe the outcome of these one-zone 
calculations (eq. 115(1: we will indicate below how our 
estimate of £ C ooi may need to be modified if sedimen- 
tation is active. We compute the ignition conditions 
as described in § 13.31 and for the first set we fix F = 
0.1 MeV m/m u . Figure shows the post-ignition evolu- 
tion for accretion rates 9.1 x 10~ 4 m,Edd (80 g cm -2 s _1 ; 
solid lines), 1.1 x 10~ 3 mEdd (100 g cm -2 s _1 ; dotted 
lines), and 2.3 x 10~ 3 m-Edd (200 g cm~ 2 s _1 ; dashed 
lines). We plot three different quantities: the evolution 
of temperature (top panel), the heat flux, F coo \ = y£ C ooi, 
normalized to the accretion flux (middle panel), and the 
mass fractions of hydrogen and helium (bottom panel). 
We then repeat this calculation at higher mass accre- 
tion rates (Fig. 0), 5.7 x 10~ 3 m E dd (500 g cm~ 2 s _1 ; 
solid lines) and 1.1 x 10 -2 mEdd (10 3 g cm~ 2 s^ 1 ; dotted 
lines). 

It is immediately clear that there are two very dif- 
ferent outcomes for the H ignition. At lower accretion 
rates, the rise in temperature following H ignition is 
sufficient to trigger a vigorous He flash with a decay 
timescale ~ 10 s (Fig. [5] middle panel). At higher ac- 
cretion rates, however, the flash is too weak to ignite 
He. As an interpretation of these results, recall that the 
ignition curve for the 3a reaction has a turning point 
at yturn ~ 5 x 10 7 g cm" 2 , which terminates the unsta - 
ble branch l|Bildstenlll99cl ICumming fc Bildstenl [20001. 
At the higher accretion rates (Fig. the ignition of H 
occurs at y < yturn, and the local rise in temperature 
does not trigger unstable He ignition. In fact, the tem- 
perature rise (see Fig. is not even sufficient to initiate 
convection, as the radiative temperature gradient needed 
to carry F cool is d\nT/d\nP w 0.25 < (<91nT/<91nP) s . 




20 40 60 80 100 
Time [s] 

Fig. 6. — One-zone burst calculation following unstable H igni- 
tion for three mass accretion rates: m/m Edd = 9.1 X 10~ 4 (solid 
lines), 1.1 X 10~ 3 (dotted lines) and 2.3 X 10 — 3 (dashed lines), 
respectively. Top panel: temperature evolution; Middle panel: the 
ratio of one-zone cooling flux to the accretion flux; Bottom panel: 
mass fraction of hydrogen (thin lines) and 4 He (thick lines), re- 
spectively. Diffusion and sedimentation is included. This category 
of H ignition triggers helium ignition and produces strong x-ray 
burst. 



In contrast, at lower accretion rates (Fig. , the ignition 
of H occurs at y > yt mn , and the rise in temperature will 
ignite the triple-a reaction. 

Thi s behavior was anticipated by iFuiimoto et alJ 
(1981), who noted that the hydrogen-shell flash would 
not develop continuously into a helium-shell flash unless 
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Fig. 7. — One-zone burst calculation following unstable H igni- 
tion for two mass accretion rates: m/mg^d = 5-7 x 10 — 3 (solid 
lines) and 0.011 (dotted lines), respectively. The temperature rise 
is insufficient to trigger He ignition 



the ac creted mass were sufficiently large. iFuiimoto et al.1 
(1981) suggested that in this case, the ignition of hy- 
drogen would lead to temperatures needed for stable H 
burning and a helium flash would be triggered, as if the 
accretion were at higher rates. We do not find this behav- 
ior. A similar case w as found in a numerical calculation 
ijAvasli &: Jossl 19*821) , in which unstable H ignition raised 
the temperature of the envelope and gradually led, after 
« 10 3 s, to unstable He ignition. This calculation used 
a much smaller neutron star radius (6.57 km). As a re- 
sult, the local accretion rate for unstable H ignition was 
much higher than ours; the reason they found an initial H 
flash was that they had assumed a very cold (5.4 x 10 7 K) 
crust and core. Unlike these previous works, we find a 
regime in which the weak H flash does not lead to sus- 
tained high envelope temperatures and a subsequent He 
flash. Because our envelope cools after H ignition, the 
accumulated He should ignite at a much greater depth 
than found previously. We discuss this further in the 
next section. 

Because sedimentation concentrates the CNO abun- 
dance into a thin layer, the expression for £ coo i, which 
implicitly assumes that the burning is over a scale height, 
may underestimate the cooling rate. We find that at 
rh = 1.1 x 1CT 3 m Ed d (100 g cm" 2 s" 1 ), the CNO is con- 
centrated into a layer of thickness Ay/y » 0.2. The value 
of £ C ooi could therefore be larger by up to (y/Ay) 2 rj 25. 
As a test, we also compute our one-zone run for this ac- 
cretion rate with y, in the expression for e coo i, multiplied 
by £ so that e CO oi is multiplied by £ 2 . We find that for 
£ < 3, there is little change in the evolution of the burst; 
at £ = 5, however, the peak temperature is reduced by 



a factor rj 2 and the burst is considerably weaker. In 
the absence of a multi-zone calculation, we cannot deter- 
mine with certainty the increase in e coo i; so long as our 
value is not more than a factor of » 10 too low, however, 
our calculations should remain valid guides to the actual 
behavior. 

We repeated the preceding one-zone calculations for 
initial conditions computed without sedimentation and 
diffusion (Figs. and EJ. The strong bursts in Fig. (HI 
have a slower rise and broader peak, but decay on a sim- 
ilar timescale, which is expected as the ignition column 
is similar in both cases. Notice that the H mass frac- 
tion decreases from r* 0.7 to rs 0.3 in the 10 4 s before 
He ignition is triggered. Although the initial instability 
is triggered by H capture onto CNO nuclei, this quickly 
saturates, and the hydrogen burning rate goes to the 
/3-limited rate ehcno = 5.8 x 10 15 Zcno ergsg" 1 s" 1 . 
The smaller CNO mass fraction (0.02 compared to 0.1 
with sedimentation included) corresponds to a smaller 
energy generation rate, which requires a longer burning 
time and more H fuel exhaustion to reach the same He 
ignition temperature. The overall amount of hydrogen 
available for the rp-process is larger, however, which ex- 
plains the broader peak. The similar decay time results 
from the similar amount of H left after the peak. For 
the same reason (smaller energy generation rate), the H 
flashes that occur at higher accretion rates (Fig. EJ) are 
much weaker than those with sedimentation and diffu- 
sion included. The hydrogen depletes on a much longer 
timescale because of the lower abundance of CNO nuclei. 
It is the abundance of CNO nuclei that sets the available 
heat deposition and hence the peak tempe rature reached 
durin g the flash, < 10 8 (Z C no/0.01) K ( Avasli & Joss 
1982). The difference between Fig. [3 and Fig. is thus 
the amount of H consumed; sedimentation enhances the 
abundance of CNO at ignition and there is a sharper 
temperature rise. 

When sedimentation is neglected, the transition accre- 
tion rate, below which the H ignition triggers vigorous 
bursts, decreases. For example, in the absence of sedi- 
mentation, the H ignition at a rate rh/rh^dd = 2.3 x 10 
does not lead to ignition of the triple-a. The reason is 
that the decreased He abundance, V( 4 He) ~ 0.7, causes 
2/tum to increase, in this case to yturn ~ 8 x 10 7 gem" 2 . 
Under these conditions, H ignites at y < ytum for rh < 
10" 3 m Ed d (cf. Fig. EH. 

We computed the ignition column for H burning for 
a range of rh and F. Figure 1101 displays the H ignition 
column density as a function of rh for Fm u /rh = 0.1 
(open squares), 0.5 (open triangles) and 1.0 MeV (open 
circles), respectively. We also plot a representative value 
of J/turn to guide the eye (horizontal solid line). A larger 
emergent flux F increases the temperature and decreases 
the H ignition column density. Notice that if we were to 
plot the H ignition column against F the curves would 
line up; the curve Fm u /rh = 1.0 MeV is approximately 
a translation of the curve Fm u /rh = 0.1 MeV shifted 
downward by 10 in accretion rate. In the absence of 
steady H burning, it is the flux emergent from the crust 
that sets the temperature structure. 



4.2.3 

In these one-zone 
the timescale for the 
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calculations of weak H flashes, 
H to be consumed is ~ 



10 



Peng et al. 








20 



40 60 
Time [s] 



80 100 



Fig. 8. — One-zone burst calculation following unstable H ig- 
nition for two mass accretion rates: m/rh^d = 9.1 X 10 — 4 (solid 
lines) and 1.1 X 10 — 3 (dotted lines), respectively. Similar to Fig.lGl 
but for cases without diffusion and sedimenation. The mass frac- 
tion of H decreases from 0.7 to 0.3 in the 10 4 s prior to He ignition. 
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Fig. 9. — One-zone burst calculation following unstable H igni- 
tion for three mass accretion rates: fa/fn-^dd — 2.3 X 10~ 3 (solid 
lines), 5.7xl0 — 3 (dotted lines) and 0.011 (dashed lines), respec- 
tively. Similar to Fig. but for cases without diffusion and sedi- 
menation. 



10 3 (0.1/Z CN o)(^h/0.3) s < y/m. We are unable to 
follow the system through a large number of repeated 
bursting cycles, but we note here that if subsequent H 
flashes do not ignite the underlying He, then a large layer 
of nearly pure He will accumulate. Because our system is 
not burning H steadily, the temperature in the He layer 
is colder than if the burning were in steady-state. As a 
result, a large He layer should accumulate. 

To illustrate this, we show in Fig.^Jthe thermal struc- 
ture of an accumulated pure He layer for rh/rfiEdd — 
5.7 x 10 -3 and 0.011. We integrated the thermal struc- 
ture equations (eq. [5] and ^Uj with d/dt — » 0) from 
the column where H ignition occurred, and set the tem- 
perature there to that found in the accumulating model. 
That is, we neglect heating of the H layer by the weak H 
flashes. In tha t case, the probl e m bec omes similar to that 
computed bv lCumming et alJ ((2,005). For our purposes, 
we treat the flux emergent from the crust as a free param- 
eter and show our solutions for Fm u /fn — 0.1 (dottted 
lines), 0.2 (dashed lines) and 1.0 McV (dot-dashed lines). 

Superposed on this figure is the ignition curve for 
pure He (solid line) defined by a one-zone stability cri- 
terion, (<9e nuc /<91nT) F = (<9e coo i/dlnT),p. The flux 
from the crust sets the temperature gradient in the ac- 
creted He, and we show the thermal structure for three 
different values, Helium ignites at y w 2 x 10 10 -3 x 
10 11 gem -2 (recurrence time « 1 - 19 yr) and at y fa 
2 x 10 9 -6 x 10 10 gem -2 (recurrence time rj 0.06 — 2 yr) 




for m/mBdd 
and 10 3 g cm 



= 5.7 x 

2 c—V 



10~ 3 and 0.011 (500 g cm-' 
, respectively. These results are 



Fig. 10. — Ignition column density of hydrogen burning as a 
function of mass accretion rates for F m u /rh = 0.1 (open squares), 
0.5 (open triangles), and 1.0 MeV (open circles), respectively. The 
solid line indicates the minimum column for the unstable ignition 
of a pure He mixture j/turn- 



sensitive to the flux fr om the deep crust, as noted by 
iCumming et all (|2~005). The flux from the deep crust 
depends on the composi tion of the ash pro ducts and 
the mass accretion rates l)Brownll200"ot |2004|) . A larger 
inner flux increases the temperature and decreases the 
He ignition column. These He flashes at great depth 
should be similar in cha racter to those of 2S 0918—549 
ijin 't Zand et alJ l2005b|) . Our ignition columns for a 
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Fig. 11. — Steady-state thermal structure of an accumulated pure 
helium layer for two mass accretion rates: rh = 500 (rh/mEdd = 
5.7 X 10~ 3 ; thick lines) and 10 3 g cm -2 s -1 (m/rh. E dd = 0.011; 
thin lines), respectively. The results for three different inner fluxes, 
F m u /rh = 0.1 (dotted lines), 0.2 (dashed lines) and 1.0 MeV (dot- 
dashed lines), are presented. The pure helium ignition curve (solid 
line) is also shown. 



given F and rh are less than that of lin 't Zand et al.l 
( 2005b) because our temperatures in the accreted H en- 
velope are higher. 

We were initially quite excited that energetic He bursts 
could also be observed from systems accreting H-rich 
material over a range of accretion rates. Indeed, we 
note that the inferred mass accretion rates listed in 
the bottom part of Table El are roughly consistent with 
the scenario we have sketched here. There is a poten- 
tial inconsistency, however. We find weak H flashes for 
m/rhEdd ~ 10~ 2 for F m u /rh = 0.1 MeV, but this low 
flux leads to huge He ignition depths (y ~ 10 11 gcm~ 2 ). 
The cooling timescale of a burst ignited at this depth 
is an order of magnitude longer than what is observed 
from the long bursts. A larger flux from the crust , which 
is mo re appropriate for low mass accretion rates l|Brownl 
l200l . increases the temperature and decreases the He 
ignition depth, as shown in Fig. El This increased flux 
decreases, however, the range of mass accretion rates 
for which weak H flashes can occur. For example, at 
Fm u /rh = 1.0 MeV, the range of mass accretion rates 
for weak H flashes decreases to 3 x 10~ 4 rriEdd rh ^ 
10~ 3 mEdd (see Fig. which is lower than the accre- 
tion rate inferred for these bursters. Given the difficulty 
in inferring mass accretion rates from observed fluxes and 
the crudeness of our one-zone calculations, it is uncertain 
whether this mismatch is a serious problem for our model 
or not. 



4.3. 



12 C Production and Implications for Superbursts 
have been detected from 10 sources 



Superbursts 
keri_ ! 2002 



to date llCornelisse et all l2000t lin 't Zand et alJ 120001 
lin 't Zand et al 



Kuu 



2001 



iStrohm aver 



Browr 
iKuulkersI I2005D . 



mean recurrence time o f these superbursts is 
ijin 't Zand et alJ l2004bj) . and the majority of super 



The 
1.5 yr 



bursters discovered to dat e accrete hydrogen-rich mat- 
ter at rates 0.1-0.3 rfiEdd- iCumming fc Bildstenl (2001) 
proposed that the unstable ignition of a small amount 
(~ 10% by mass) of 12 C at densities < 10 9 gem" 3 
would suffice to power these bursts. This scenario pre- 
dicts recurrence times, energe tics, and cooling timescales 
ijCumming fc Macbeth! |2004|) that are in good agree- 
ment with observations provided that t he deep crust 
of the neutron star is sufficiently hot jBrownl 120041: 
iCooper fc Naravanl2005HCumming et al.ll2005j) . A chal- 
lenge to this scenario is the lack of 12 C produced 
in both one-zone and multi-zone ca lculations of X-ray 
bursts with full reaction networks iSchatz et alJl200H 
i Wooslev et all 120041: iFisker et all l2005al iFisker et alJ 
2005b), although a greater abundance of carbon may be 
produc ed if the burning at lower accretio n rates were 
stable ljCumminell2005l: ICooper et alJl2005|) . 

As shown by Fig. sedimentation decreases the H 
abundance and increases the CNO abundace at the 
base of the accreted envelope, even at accretion rates 
> 0.1 rh-Edd- It is the ratio of protons to heavier "seed" 
nuclei that determines th e run of the rp-pr ocess (for 
a succinct discussion, see iSchatz et all Il999j) . In the 
absence of mixing, sedimentation lowers the proton-to- 
seed ratio and decreases the mean mass of the rp-process 
ashes. As we noted in § 14.21 during the burst rise a 
convective region is established which moves outward 
in mass and mixes the envelope; this prevents us from 
giving a definitive statement about how the stratifica- 
tion of the envelope changes the burst ashes. Just to 
motivate the problem, we perform a one-zone calcula- 
tion (§ 14.21 eq. 22]) of the burst nucleosynthesis. The 
initial temperature, column density, and composition 
are taken from the values of the quasi-static calcula- 
tion at the ignition point for a mass accretion rate of 
rh = 0.11 TOEdH- We checked our calculation against the 
results of Sch atz et alJ i|2001fk for similar conditions, we 
find a broad distribution of isotopes with A — 60-100, 
but our distribution peaks at A ks 64. 

Figs . 1 1 21 and 1131 show the results of this one-zone calcu- 
lation, for ignition conditions calculated both with (solid 
lines) and without (dotted lines) a depressed hydrogen 
abundance. In Fig. El we plot e coo i (a crude measure of 
the burst lightcurve; top panel) and the mass fraction of 
12 C, while Fig. 1131 shows the composition of the ashes. 
The temperature rises rapidly for t < 5 s during the he- 
lium flash. Helium burns via triple- alpha reaction and 
produces car bon, which i s cons umed via 12 C(p, 7) 13 N. 
As noted bv ISchatz et alJ (|2003), 12 C is only abundant 
after the H has been consumed. When X(H) is decreased 
by a factor w 2 (Fig. EJ, 12 C is produced earlier and 
in greater abundance. In addition, the lower proton-to- 
seed ratio shortens the burst and shifts the composition 
to lower mass numbers. For this accretion rate, A( 12 C) 
is increased by a factor of 5 (0.01 to 0.05) in the ashes. 
We repeated this calculation for rh — 0.057 rfiEdd and 
rh = 0.23 rriEdd, and found that A( 12 C) increased by a 
factor w 3 and ks 7, respectively: the overproduction of 
12 C is greater at higher accretion rates. 

As a check, we also performed the run at rh = 
0.11 mEdd with varying resistance coefficients K^. As 
disc ussed in Appendix[A] the Kij are very uncertain (see 
Fig. IA14|) . Reducing by a factor of 5 enhances the 
sedimentation of He and CNO elements and leads to an 



12 



Peng et al. 



increase in the 12 C mass fraction after the He flash by 
a factor of ~ 2 (X( 12 C) = 0.09). Conversely, increasing 
Kij by a factor of 5 (smaller Dij), decreases X( 12 C) by 
a factor of 2 (X( 12 C) = 0.03). 




10 20 30 
Time [s] 



50 



Fig. 12. — The one-zone cooling rates e coo l (cq. 1151 1 and the mass 
fractions ol 12 C for a one-zone burst calculation at m = OAlrh^dd- 
The initial composition is taken from the bottom of the fuel layer at 
the point 4 He ignites, for the calculation with (solid) and without 
(dotted) sedimentation and diffusion. 




40 60 
Mass Number 

Fig. 13. — The ashes produced by the one-zone burst calcula- 
tion following unstable He igntion for m = GAlm^d- The initial 
composition is taken from the bottom of the fuel layer at the point 
4 He ignites, for the calculation with (dots, solid lines) and without 
(crosses, dotted lines) sedimentation and diffusion. 



5. SUMMARY AND CONCLUSIONS 

Using a simplified numerical model of an accreting neu- 
tron star envelope that allows for differential isotopic ve- 
locities, we have undertaken a first study of the effects 
of sedimentation and diffusion on the unstable ignition 
of hydrogen and helium on the surfaces of accreting neu- 
tron stars, and have investigated the outcome of unsta- 



ble hydrogen ignition using simple one-zone models. Our 
principal conclusions are: 

1 . The effect of sedimentation changes the conditions 
at H/He ignition even for to > O.lmEdd- 

2. There is a range of accretion rates for which un- 
stable H ignition does not trigger a He flash. This 
range depends on the flux from the deep crustal 
heating and the degree of settling of He and CNO 
nuclei. For F (to u /to) = 0.1(1.0) MeV, the range is 
3 x 10~ 3 m Ed d < rh < 10~ 2 m Edd (3 x 10~ 4 TO Ed d < 
to < 1 0~ 3 TOF.dd). In contrast to previous calcu - 
lations l|Fuiimoto et alJll981t lAvasli fc .Tosslll982lK 
we find that the flash does not lead to quasi-steady 
H burning that develops after « 1000 s into a He 
flash. We speculate that successive weak H flashes 
can lead to the accumulation of a large reservoir of 
He; this may explain the long bursts observed from 
some sources. For accretion rates lower than this 
range, H ignition leads to a strong mixed H/Hc 
flash. The rise time of the flash and the burst 
duration, however, depend on the abundance of 
H and CNO nuclei at the ignition. The sense of 
this (a transition from vigorous bursting to weak 
H flashes interrupted by less frequent, strong He 
flashes) agrees with the observations of bursts at 
low m, although our estimates of the accretion rate 
at which this transition happens does not match 
the observed value. 

3. For to > 0.1 TOEdd, the partial stratification of the 
envelope means that the base of the accreted layer 
is deficient in hydrogen, relative to the case where 
complete mixing is assumed. Because a convective 
zone is established during the onset of the burst, 
we cannot say for certain what effect, if any, this 
will have on the outcome of the burst. We do note, 
however, that in the absence of mixing, the pro- 
duction of 12 C is enhanced, by a factor of 5. 

There are several items that are open for further 
investigation. First is the interplay between convec- 
tion and the rp- process breakout. At these accretion 
rates, the unstable burning of 4 He and the breakout 
from the HCNO cycle drives a convective zone out- 
wards that encompass es > 90% of the accreted envelope 
l|Weinberg et aUl200fl) . This would rapidly mix the en- 
velope and "reset" the amount of protons available for 
capture onto the endpoint of the (a,p)(p, 7). Observa- 
tionally, the onset of stability is closer to to w 0.3 TOEdd 
and not TOEdd as found by most one-dimensional cal- 
culations. In superbursting sources, much of the ac- 
creted fuel must be consumed stably to account for the 
low values of a = (per sistent fluence) / (burst fluence) 
ijin 't Zand et alJ l2004tfl . Close to stability, the 
growth of the convection zone is reduced (see 
iHeger. Cumming. fc Wooslevl I2005L for a discussion of 
the behavior near stability; note that in their calculation 
this occurs at close to the Eddington limit). The extent 
to which sedimentation affects the nucleosynthesis in this 
regime therefore remains an open question. A second 
task is to perform one-dimensional calculations of unsta- 
ble H ignition under conditions for which the ignition 
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does not trigger a He flash. This is critical for determin- 
ing how sedimentation modifies the finite-difference ap- 
proximation for the cooling rate, as described in § 14.2.21 
Following successive weak flashes is necessary to deter- 
mine how large a mass of He can in fact be accumulated. 
Finally, as discussed in § 14.31 our results depend on the 
values of resistance coefficients. A further investigation 
of resistance coefficients at T > 1 will help refining our 
results. 
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APPENDIX 

RESISTANCE COEFFICIENTS 

In this appendix, we discuss our choice of the resistance coefficients Kij for a multicomponent plasma. When the 
plasma is sufficiently rarefied, particle pairs interact via a screened Coulomb potential, and the resistance coefficient 
in this limit is 

23/2 . . / ^ \ 1/2 



K° 



-HiUj ln(l + ITjA 



Here /iy = AiAj / (A; + Aj) and Xi 



fc B T/(47re 2 Zfm)] 1/2 being the Debye length. 



3 " r ' 13 \k B T J k B T 

A\ Ji k B T/{e 2 Z l Z ] ), with A D = [ 
The sum in the expression of Ad is over all species in the plasma, including electrons. Using 



(Al) 



r = 



ak B T 



(A2) 



and the ion plasma frequency 



we rewrite equation (|A1|) in terms of T and 



(A)m n 



K 



(67T) 



-1/2 



ZiZj \ 

Here A = ln(l + a;? ). This form of K^ is for our definition for T 



(A) 2 
AjA-, 



(A) 



1/2 



^ P pr 3 / 2 A. 



(A3) 



(A4) 



(eq- 



When the plasma becomes strongly coupled, Ad < a or equivalently, when Y > 1, the resistance coefficient can no 
lon ger be defined in terms of a sequence of binary collisions. This regime is important for studies of white dwarf cooling 
fsee lDelove fc Bildstenll2"002l) , and there have been numerous att empts to calculate the bin ary diffusion coefficient Dij 
from molecular dynamics simulations (for a recent example, see iDalig ault fe Murillo|| 2005j) . The resistance coefficient 
is related to by — (k B T '/ ' Kij){riinj /n), where n — Y^ s n.. lHansen et alJ |)198It(1 found that 



Kij = o.iipwpr' 



0.34 



(A5) 



iBildsten fc Halll II2001D used the S tokes-Einstein relation with the shear viscosity estimated from molecular dynamics 
simulations ijDonko fc Nvirill2000() to obtain a similar relation, 

K i:j = 0.15pcj p r - 3 . (A6) 

Both of these fits are applicable in the regi me T > 1. For a m ulti-species plasma, there are no well-defined binary 
diffusion coefficients in the regime of T > 1. lHansen et alJ l)1985fl noted that a linear combination of the self-diffusion 
coefficients computed for a OCP gave a first-order accurate estimation of the binary diffusion coefficient Dy; when 
written in terms of the resistance coefficients, this becomes 



Ki 



K K 



(A7) 



In applying this formula, we also use a p in place of a, where Aira^/i = (Z)/n e in computing the Ku. This follows 
from the Stokes-Einstein relation ijBildsten fc Hal]|l2001() . 

The long-range cutoff for the Coulomb potential in equation QA4|I does not hol d when a > Ad- iMuchmorel (|l984) 
proposed redefining x%j to use max(Ao, a)', using this definition and equation (|A2|> . we have 



(z 2 ) i 
'ZiZj r 

I (z 2 ) 1 
l ZiZj r' 



mr+W>7WT> 



1/2 



Ad > a 
Ad < a 



(A8) 
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Fig. A14. — Normalized resistance coefficients for a H/He mixture with X ^H) = 0.7 and X( 4 He) = 0.3. We show the result for the 
low-density limit Kfj (dotted line) as well as the calculation of lMuchmord 119841) (dashed line; M84), the fit of llben &: MacDonald 11983) 
(solid line; IM85), and the Stokes-Einstein relation with the viscosity inferred from OCP molecular dynamics simulations IBildsten fc Hall 
2001, dot-dashed line; BH01). 



The two expressions match at T = \(1 + (Z) / (Z 2 )) -1 . The fa ctor 1 + (Z)/{Z 2 ) accounts for Ad including the electrons 
in the sum over charged species. iFontaine fc Michaudl l)1979f) numerically computed binary diffusion coefficients using 
a Thomas-Fer mi potential at high densiti es Sy < 1 and a Debye-Huckel potential at low densities. A convenient fit to 
their results is l|Iben fc MacDonaldllToS^) 

In [1 + 0.187694^1 n 
K = 1.6249 1 ^ K o Ag 

In (1 + 4) 

where computed from equation (|A8|) . 

Figure IA14I shows the normalized resistance coefficients -Ky (w p /9) _1 as a function of V for a two- isotope plasma 
with X(H) = 0.7, X( 4 He) = 0.3. We show Kij computed from a Debye-like potential with a screening length 
max( Ar), a) ijMuchmorel 119841 dashed line); from a Thomas-Fermi potential at higher densities l)Iben fc MacDonaldl 
1985, solid line); and fr om the Stokes-Einstein relation with the viscosity inferred from molecular dynamics simulations 
(Bild sTen fc Hallll2001l dot-dashed line). We also show (dotted line) for comparison. In the regime 0.1 < T < 10, 
the fit of llben fc MacDonaldl l)1985j) agrees well with the calculation of lMuchmord ll 19841). and has the correc t scaling 
with r in the low- and high-dens ity limits, where it diffe rs by < 40% from the fit of IBildsten fc Hafil l|2001f ). It also 
agrees well with the calculation of lPaauette et all l)1986al) over the range 0.1 < T < 10.0. For the conditions of interest 
in this paper, generally T < 1, we therefore adopt this fit throughout this paper, with the caveat that the resistance 
coefficients are generally uncertain in the range T > 1. 
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